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We present in a simple terms the theory of the factorizing operator introduced re- 



cently by Maillet and Sanches de Santos for the spin- 1/2 chains. We obtain the explicit 
expressions for the matrix elements of the factorizing operator in terms of the elements 
of the Monodromy matrix. We use this results to derive the expression for the general 
scalar product for the quantum spin chain. We comment on the previous determination 
of the scalar product of Bethe eigenstate with an arbitrary dual state. We also establish 
£N) ■ the direct correspondence between the calculations of scalar products in the F - basis 
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and the usual basis. 



1. Introduction. 



43 

One of the most important open problems in the theory of quantum integrable models is 
the calculation of the correlation functions. In the framework of the Algebraic Bethe Ansatz 
method |J the problem is the combinatorial complexity of calculations due to the structure 
of Bethe eigenstates. Only the very limited number of physical results for the correlation 
functions (or formfactors) of integrable models have been obtained from the first principles. 

The concept of factorizing F - matrix was recently introduced by Maillet and Sanches 
de Santos M following the concept of Drinfeld's twists in his theory of Quantum Groups. 
The existence of the factorizing matrix allows one by means of similarity transformation to 
make the elements of the Monodromy matrix totally symmetric with respect to an arbitrary 
permutation of indices 1,...N together with the corresponding inhomogeneity parameters 
£1, . . . £/v- The matrix Fi_n is denned as follows. Consider any permutation a G Sjy Then 
any component in the auxiliary space of the monodromy matrix T (t) depending on the 
parameter t transforms as 

To >(T i...(tN — Ri...nTo,i...n{Ri.,.n) 1 '> 
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where R° N is the operator constructed from the elementary S- matrices. Then if the matrix 
Fi . jv is defined according to 

(-^l...Af) _1 = (Fl...N)~ 1 Fal...aN, (1) 

in the F- basis the monodromy matrix Tq = FT F~ l is totally symmetric: 

T^al-aN^) = ^l...jv(*)> CT E S N 

for arbitrary t. It was also realized that there exist the factorizing matrix F which diagonalize 
the A (D)- component of the monodromy matrix. In this basis the operators B and C have 
quite a simple quasilocal form which recently made possible the direct computations for the 
important class of the (general) correlation functions 0. In some respect this approach allows 
one to simplify the calculations in comparison with the general theory of scalar products [|J, 
|| and gives an alternative way to obtain many important results for the XXZ spin chain. 
Among them one can mention the derivation by Korepin the Gaudin formula for the norm 
of Bethe eigenstates Q and the calculation of the scalar product of Bethe eigenstate with an 
arbitrary dual state |6[], which leads straightforwardly to the determinant representation of 
the formfactors of basic (local) operators (see e.g. 0, || for the analogous calculations for 
the case of Bose-gas with S - function interaction). Using the Algegraic Bethe Ansatz and the 
solution of the quantum inverse scattering problem the authors of ref . || obtained the multiple 
integral representation for the correlation functions found previously using the other methods 



(see for example |L0| and references therein). 

At the same time in the approach of ref.0 the construction of the F - matrix itself and 
of the matrix elements of different operators in the F - basis is quite complicated. The 
construction involves the specially defined partial F - matrices and some steps (for example 
the proof that the factorizing matrix diagonalize the A (D)- operators) are based on a special 
recurrence procedures. Also the explicit expressions for the matrix elements of the operator 
F, which can be used in practical computations, was not presented. The generalization of the 
construction to the case of the other models is not straightforward. 

The aim of the present paper is twofold. First, we present in a simple terms the theory 
of the factorizing operator introduced in ref. and derive the new formulas for the matrix 
elements of the factorizing operator (Section 2). In comparison to the authors we proceed 
in a different way. First of all we construct the matrix that diagonalizes the operator A(t). We 
obtain the simple expression for its matrix elements and the matrix elements of the inverse 
operator in terms of the matrix elements of the products of the operators B(t), C(t). These 
expressions have not been obtained previously. Then we find the expressions for the other 
operators B, C in the new basis (Section 3). Afterwards we prove that the constructed 
operator is the factorizing operator as defined above. This procedure allows one to obtain 
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the simple expressions for the matrix elements of the factorizing operator which can be used 
in practical computations. The presented formalism could be usefull for the other integrable 
models (for example the generalization to the case of XYZ- spin chain is straightforward). 

Second, we apply the developed formalism to the calculation of the general scalar product 
in a way which is different from that of ref. || (Section 4). We obtain the general expression 
for the scalar product of ref.pl] in a direct (i.e. without any recurrence procedures) and a 
simple way. In some sense our calculation clarifies the mathematical structure underlying 
the derivation of ref's ||, @- The key point in our derivation is the formulas for the matrix 
elements of the factorizing operator obtained in Section 2. Next, in Section 4, we supplement 
the derivation || of the determinant representation of the scalar product of Bethe eigenstate 
with an arbitrary dual state. We present here the new version of the proof given by Slavnov 
0, which is based on the recurrence relations. We also establish the direct correspondence 
between the calculations of an arbitrary scalar products in the F - basis and the conventional 
basis in the Appendix C. The calculation of the sums for the scalar products presented here 
can serve as an independent proof of the formulas for the matrix elements of the Monodromy 
matrix in the F - basis. 

2. Construction of the factorizing operator. 

We consider in this paper the XXX or XXZ spin- 1/2 chains of finite length N. Before 
dioganalizing the operator A, let us fix the notations: the normalization of basic S - matrix, 
the definition of monodromy matrix and write down the Bethe Ansatz equations. For the 
rational case (XXX- chain) the S- matrix has the form S^^i, £2) = ti — t 2 + tjP^, where P\ 2 
is the permutation operator. In general (XXZ) case it can be written as 



Sl2(£l, £2) — 





a(t) 








\ 







c(t) 


m 










bit) 


c(t) 





V 











a(t) ) 



One can choose the normalization a(t) = 1 so that the functions b(t) and c(t) become: 

b(t) = 4>(v)/4>(t + v) c(t) = 0(t)/0(t + v ), 

where <ft(t) = t for the isotropic (XXX) chain and </>(£) = sinh(t) for the XXZ- chain. With 
this normalization the S'-matrix satisfies the unitarity condition S\2{t\, t2)5 , 2i(^2, £1) = 1- The 
monodromy matrix is defined as 

T (t,{Q) = S 10 (^,t)S 20 (^,t)...S N0 ^ N ,t), 
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where & are the inhomogeneity parameters. We define the operator entries in the auxiliary 
space (0) as follows: 

\To\a) = I *® ~S I : "• ? = (1.2) = (|:j). 



C(f) £>(t) 



a/3 



We denote throughout the paper (T; !) = (1 ; 0) so that the pseudovacuum (quantum reference 
state) |0) = |{00...0}jv). The triangle relation (Yang-Baxter equation) reads: 

S12S13S23 = S23S13S12, RqO'TqTqi = TqiTqRqqi] Rqqi = Poo'Sqo'- 

The action of the operators on the pseudovacuum is: A(t)\0) = a(t)\0) (a(t) = n a c(£, a — £)), 
D(t)\0) = |0), C(t)\0) = 0. The Bethe Ansatz equaions for the eigenstate of the Hamiltonian 
Yltx B(ti)\Q) are 

and the corresponding eigenvalue of the transfer - matrix Z(t) = A(t) + D(t) is 

M M 

A(t, {t a }) = a(t) J] c-\t a - t) + J] c'\t - t a ), 

a=l a=l 

where t a are the solution of Bethe Ansatz equations. 

In order to construct the operator O = Oi...jv which diagonalizes the operator A(t) 
(0~ l AO = diag(A)), let us first construct the eigenfunctions of the operator A. One can 
do it in two different ways. First, note that A - is a triangular matrix in a sense that it 
makes particles (the spin-up - coordinates) move to the right on the lattice 1...N. Thus its 
eigenvalues are coincide with its diagonal matrix elements and therefore are charactarized by 
the set of integers (rii, ■■■Um) - the spin-up positions. Let us denote this eigenfunctions by 
\<j>(ni, ...tim))- Clearly, |{00..0}jv_m{11..1}m) is an eigenstate of A(t). Therefore, considering 
the permutation 

Sio...5jVQ — > Sio...SNoS nM o...S ni o (rii < ri2 < ... < TIm), 

we realize that 

\(f)(n 1 ,...n M )} = T ni T n2 ...T nM \n 1 ,n 2 , ...n M ) (2) 

where 

T n = S n +i^ n S n +2,n---SNn-> 

is an eigenstate of the operator A(t). Note that if we modify the given permutation inter- 
changing rii ^ n j the operator in the right-hand side of (Q) modifies but the state 
remains the same. For example, for M = 2: T n2 T^ t = T ni T n2 S~ 2 l nV where ri\ < n 2 and the 
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prime means the absence of the term S nin2 in T ni . Since S nin2 \ni,n2) = \nx,n2) the state 
remains the same. 

The second way - is to consider the state 

|0(th, ...n M )) = B(U)B(U)-B(Hn,M> (rii ? n 3 ). (3) 
Using the fundamental commutation relation: 

A(t)B(q) = -±— B {q)A(t) - ^L^B(t)A(q), (4) 

and the fact that for the pseudovacuum state A(^)|0) = 0, we find again that \4>{{n})) 
({n} = {ni, ...tim}) is an eigenstate of A(t) with the eigenvalue A{ n y{ n y(t) = I\ a ^n k c(£a — t), 
and the operator A(t) in the new basis has the following form: 




(5) 



c(£ a -t), a^n k 
1, a = n k 

One can see that the states |0({n})) (||) and @ are coincide. This can be seen using the 
identity 

B(Q=T n (0\S l0 ...S n ^,oP n o\l) (6) 
and ordering the (commuting) operators -B(£ n J in eq.@ according to the prescription rii < 
n 2 < ... < n M) so that the second operator in the last formula simply creates the particle 
(spin- up) at the site n with the amplitude equal to unity. 

Now we can introduce the diagonalizing operator O = 0\,,,n (we will see later that it is 
also the factorizing operator). Let us define the operator O such that 

\<P{{n})) = 6{{n})\{n})=6\{n}), 

where 0({n}) is given by the equation (0). Clearly the operator A F (t) = 0~ 1 A(t)0 - is 
diagonal (see @). Now we immediately get the simple formula for the matrix elements of the 
factorizing operator: 

6 {m}{n} = ({m}m ni )B(U)...B(Cn M )\0), (7) 
where {m} = {m 1; ...trm }• From this expression (see (||)) one see that O is the triangular 
matrix (in the same sense as the operator A((t)). Let us show that the operator O is invertible 
and construct the inverse operator. To find this operator we have to consider the following 
dual states: 

(4>(n u ...n M )\ = (0\C(U)C(U)-C(£ nu ). (8) 

Define the operator O analogously to the previous case as (4>{{n})\ = ({n}\0. The matrix 
elements of this operator are equal to 

{m}{n} = (0\C(U)C(U 2 )-C(U M )\{n}). (9) 
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Let us calculate the following scalar product 

(0({m})|0({n})) = ({m}|O6|M)= (10) 

(o\c(u)C(U)---C(UJB(c ni )B(u)--MUM- 

This scalar product can be calculated using another well known relation relation between the 
elements of the monodromy matrix following from the Yang-Baxter equation: 

[B(q), C(t)] = (D(q)A(t) - D(t)A(q)) , (11) 

and using again that for any site i, A(^)|0) = 0. Moving the operators A and D to the right, 
repeating consequently the relation (ffl) and the analogous relation for D (which differs from 
eq.(|J) by the interchange t <-> q in the coefficients between the products of the operators) and 
using the equations 

C(£ n )B(£ n )\0) = (n^a-U) |0>, 

we obtain that the matrix / is diagonal : 

<0({m})|0({n})> = ({m}\d6\{n}} = <{ro}|/|{n}> = (j[5 ni ^ f({n}), 

The corresponding diagonal matrix elements can be found either using the procedure men- 
tioned above or, which is the simplest way, using the representation (§) (and the same for the 
operator C(t)). In fact, taking into account the triangularity of O and O we obtain: 

({ n }\f\{ n }) — ^2 d{n}{m}0{m}{n} = 0{ n }{ n }6{ n }{ n } . 
{m} 

The diagonal matrix elements can be easily calculated using the representation (|7|) and the 
formula @. We get: 

f(n 1 ,...n M ) = Y[( I] (12) 
Thus we obtained the inverse matrix O^ 1 : 

do = /, 6- 1 = r x b. 

Before proceeding with evaluation of the matrix elements of the other operators in the new 
basis, let us mention some usefull properties of the operator O, and prove that that, in fact, 
it is the factorizing operator in a sense of the defenition (0): O = First, O and O" 1 

are the triangular matrices (upper triangular as A(t)). Second, the pseudovacuum state is an 
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eigenstate of O {O l ) with the eigenvalue equal to unity. In general we have the following 
equations for arbitrary number of particles n: 

d|{oo..ov_ n {n..i} n ) = 1 Koo.-OV-nln.-iU 

({ll..l} n {00..0} J y_ n |d=l ({ll..l} n {00..0} J y_ n | J 

and the same formulas for the inverse operators. From this formulas one can already suspect 
that O is the factorizing operator. Indeed for the particular permutation (|2|) er({n}) the 
factorizing condition is represented as 0(O a ^ n ^)~ 1 = 0({n}), where 0({n}) = T ni ..T nM and 
is fulfilled at least for the state \{n}) since due to the last formulas O tT ^)- 1 \{n}) = \{n}). 
The rigorous proof goes as follows. We construct the operator that acting on the state |{n}) 
produces the state 0({n})\{n}) . It is easy to see that the operator that fulfills the above 
requirement is: 

6 = F l F 2 ...F N , Fi = (l-fiij+Tihi, (13) 

where hi is the operator of the number of particles (spin up) at the given site i. The operators 
Fi entering ( Jl3| ) do not commute and their ordering in eq.flTB]) is important. To prove the 
factorizing property of this operator it is sufficient to consider only one particular permutation, 
say the the permutation + 1) since all the others can be obtained as a superposition of 
these ones for different %. We will show in the Appendix A that 

6 = Si+md^v, (14) 

Evidently, in contrast to (|l|), for any transmutation a G Sn we will obtain only one operator 
R\ N on the left of the operator O. Thus it is proved that O is the factorizing operator in a 
sense of eq. ([!]). 

3. Construction of the matrix elements. 

In this section we calculate the matrix elements of B(t) and C(t) - operators in the F 
- basis: B F (t) = FB(t)F~ l = 0^ l B{t)0 (and the same for C(t)). The general scheme to 
perform the calculations is to use the formalism developed in the previous section, which leads 
to the following chain of equations: 

B(tM{n})) = B(t)6\{n}} = OB F (t)\{n}} =^<t>(z,t, {n})6\{n} } x}. 

X 

where \{n}, x) is the state corresponding to the new set of coordinates with an extra spin-up 
at the site x. Thus acting by the operator B(t) on an eigenstate (0) 

B(t)B(U)-B(UM =J2(B(QB(£ ni )...B(!i nu )\0))cP(x,t,{n}), (15) 

X 
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we see that <p(x, t, {n})- is exactly the matrix element in the new basis. To get the single term 
in the sum in eq . fll5|) , we act by the operator A(£ x ) (x ^ n k ) at both sides of this equation. 



Using again the property A(£j)|0) = and eq.([|) we get for the left-hand side of ( |T5| ) 
A({; x )B(t)B(U)...B(Zn M )\0) = - (b(t - Q/c(t - £,)) B{Z x )A{t)B(Z ni )...B(£ nu )\Q), 



while for the right-hand side we get the single term with B(£ x ), which can be easily evaluated 
using again eq.(£§) and the formula 

A(QB(Q\0) = f n 5(6* - £■) ) B(Q\0), 

\a^x J 

which can be proved by direct computations. After the cancellation of similar terms at both 
sides of eq. ( |1"5| ) we get 



[x 



t, {n}) = - (b(t - Q/c(t - £,)) J] ^ -t)(l[ c(£ a - q) 



Then, using the equality 

C[t - Cx) 

we finally obtain in the operator form: 

s F (t) = 5>j m -t)u( 5(a " tm « ~ ^ Q ^ n * V (i 6 ) 

a ^x V 1. a = n k J 

With this expression one can prove the equation Yli B F (£ ni )\0) = \{n}) which is consistent 
with the formulas of the previous section. 

For the operator C(t) proceeding in a similar way and using the relation ([0]), we get 

C(t)B(U)...B(tn M )\0) = B(U)C(t)B(U)-B(tn M M 

- (Kt ~ U)/c(t - U)) D(U)Mt)B(U)-B(UM- (17) 

Using the symmetry between m, ■■■Um-, we can concentrate on the term which describes the 
flipping of the spin on the n\- site and consider only the second term in the right-hand side of 
the equation (fPTD. Commuting A and D operators to the right in (|I7D and denoting £ m = £ a 
we obtain similarly to the previous case 

c F (t) = e °x - 1) n [ ~f ( : ~ t] : , w a ^ nk ). m 
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The operators (|T6|) and ( |18D are quasilocal i.e. they describe the flipping of the spin on a 
single site with the amplitude depending on the positions of spins on all the other sites of 
the chain. The operator D F (t) can be found using either the same method or the quantum 
determinant relation and has a (quasi)bilocal form. 

4. Calculation of the scalar products. 

In this section we use the developed formalism to obtain the expressions for the correlation 
functions for the spin chains. We use the factorizing operator and, in particular, the expression 
for its matrix elements (0) to obtain the expression for the general correlation function (scalar 
product) 

S M ({X},{t}) = (0\C(X l )C(X 2 )...C(X M )B(t 1 )B(t 2 )...B(t M )\0) 1 

where {A} and {t} are two arbitrary sets of parameters (not necessarily satisfying the Bethe 
Ansatz equations). The correlation function can be represented in the following form 

(0|Tr 0l ,...o 2M (°o 1 -- (X L a o M+l -- (J o 2M T Oi ■ ■ -To 2M ) |0) 

The auxiliary spaces Oi, . . . 2 m can be considered as a lattice consisting of 2M sites with the 
corresponding spectral parameters Ai, . . . Am, t±, , . . . £m- Rearranging the basic S- matrices 
entering the product of the monodromy matrices we arrive at the operator 

T\. . . T/v, 

where the new monodromy matrices act in the auxialiary space instead of the original quantum 
space: 

T n = S n o 1 S n o 2 • • • S n Q 2M . 

Obviously using this matrices the correlator can be represented as the following matrix element 
in the new quantum space Oi, ...0 2 m'- 

({00..0} A/ {ll..l} A/ |i(6) • • • i(6v)|{ll..l}Af{00..0} A/ > 

(we use the symmetry A(t) <-> D(t), <->• 1 here). Transforming the operators to the F 

- basis we find at & = 0: 

Sm = E({00..0} m {11..1}m|O|M)(M| (^(0)))" |{n})({n}|6- 1 |{ll..l} A/ {00..0} M ). 

M 

The sum is over the states labeled by the positions of M particles {n} = ni...nM on a 
lattice consisting of 2M sites with the inhomogeneity parameters Ai . . . Xm, h . . . t\t- We use 
O^ 1 = f~ l O and the representations (|7|), (Jty for the matrix elements in the last formula. Let 
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us denote by fi%, . . . hm the parameters (from the set ({A}, {t})) corresponding to the sites 
rii . . . tim, and by u±, . . . vm the rest of the parameters so that {A} U {t} = {fi} U {u}. The 
first matrix element in the formula for Sm is equal to 

({00..0} M {ll..l}A/|5(/i 1 ) . ..B{ m )\0) = ({U..l} M \B'(p 1 ) . . . B'{ m )\0), (19) 

where B'(^i) are the same operators difined on the lattice consisting of M sites with the 
inhomogeneity parameters t\ . . . tu- The second matrix element can be reduced to the same 
expression with the parameters Ai . . . Am (using the symmetry C «-> B, 1). Then using 
the formula for the matrix elements of the operator f^ 1 (see eq. flT^) ) we finally obtain the 
expression: 

s M ({x},{t})= e ( n «"s) ) a n ~ ( „ \ „ v ( 2 °) 

ni,..n M \ j / i,j V i> 

where we denoted by $m(£, t) the functions in the right hand side of eq.(ID). The determinant 



representation of this function is well known (see for example H, |]TT 



^ *) = n ^"nn^i ^ det ( ft TW p (21) 

for the rational case. Here {£} are the inhomogeneity parameters and {t} are the arguments 
of B- operators (see eq.([Lj])). We present the basic properties of this function used in the 



calculations in the Appendix B. The functions a(y) in eq.(|20"D are exactly the functions defined 
above a(v) = Yl a c(£ a — v) while in the rest of this formula due to the definition of the matrices 
Tj one should interchange the arguments in the functions c~ x {yi — fij) which is taken into 
account in eq. (|2Cf ) (or one could make the replacement i] — > —77). Using the properties of this 
functions one can represent the general formula ( p0|) in a different way: 

M 1111 



EE(II (a(A„)a(4))) $ m (t fc , \f>)& M -m(K, t a ) - 



m=0 k,n 



C(\f3 ~ A n ) c(t a - t k ) C(t a - A n ) C(\f3 - t k ) 



where the sum is over the two sets ki, ...k m and n 1; ...n M _ m . We used the following simpli- 
fied notations in this formula. We devided the set {t} into two subsets {t} = {t k } U {t a } 
where {t k } = (t kl ...t km ) G {v}, {t a } G {/j} and analogously {A} = {A„,} U {A^}, {A n } = 
(A ni ...A„ M _ m ) G {v}, {Xp} G {fi}. The products in the last formula are over the indices la- 
beling the elements of the corresponing sets. One can further rewrite this formula using the 
expressions for the functions $ m to get the formula which leads to the determinant represen- 
tation for Sm after the special dual fields are introduced ||. 

Let us derive the scalar product of the Bethe eigenstate with an arbitrary dual state along 
the lines of ref.|| starting from the formula (P0|). The set of the parameters {t} obey the 
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Bethe equations so for each a(ti) in the sum (|20|) one can substitute the function 

f (+\ = TT ^ ta - ^ = TT -tj-y) 

Then the idea is that one can calculate the sum in (p0|) for an arbitrary smooth function a(A), 
which behaves at least as a constant at infinity, used for the terms a(Xj) (not necessarily equal 
to a(A) = Ila c(^ Q — A)). In that case the sum has the simple poles in the variables {t} and 
{A} at the points £j = Xj and don't have any other poles, for example the poles at Aj = Xj, 
which exist if, according to || one considers a(Xj) as an independent variables. These poles 
was not considered in || (actually only the behaviour of the function in the parameters tj 
was considered - obviously there is no poles at £, = tj if the Bethe Ansatz equations are taken 
into account). One can use the symmetry of Sm in {t} and {A} and single out the residue at 
Ai — > ti (which is contained in the term with Ai G {u}, t\ G {/i} and vice versa): 

S M ({A}, {*}, afa)) U^ fl - ^-[^ n w . 1 . 1 rr ^M-i ({A}', {*}', a'fo)) , 

where {A}', {£}' do not contain Ai, ti and 

// \ / \C(^j — ti) 



c(*i - I/j-)' 

The functions /(tj) entering the sum Sm-i are also modified so that the terms with t\ are 
absent in their definition. Also the behaviour of Sm ~ 1/Ai at Ai going to infinity should be 
taken into account. It can be easily proved that the same recurrence relation at Ai — > t\ is 
obeyed by the following expression: 

S M ({X},{t}) = n t, ,L ft rT det(M,,(t,A)) 

\\i<j\ti - tj) Uj<i{Xi - Xj) n 

Mij(t, A) = _\ I a(A,) I](*a - A, + V ) - H(t a - X, - V ) ) , (22) 

where at the end of calculations one should take the actual form of the function a(A). So we 
proved that for an arbitrary function o(A) the residues of the only existing poles at U = Xj 
satisfy the the same recursion relation in both formulas. Since the function (in each variable) 
is completely determined by the positions of (all) poles and the corresponding residues (along 
with the behaviour at the infinity) the two functions fl2"2"|) and fl2"0|) should coincide for an 
arbitrary function a(A). To complete the proof, one should check the formula (|22| ) for M = 
1,2. Equivalently, one can check the coefficients corresponding to the terms Il^i a (Ai) and 
ni=i a (U) = 1 in eq.(p2]) which are respectively 
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in agreement with ( pO|) . The formula (22) is written for the rational case. The generalization to 



the case of XXZ chain is obvious. This formula can also be represented through the Jacobian 



as: 



. i/ Ui,j(U - Xj) .!.„ ( <) 



S M ({X},{t}) = n n Yn n det \W A ^ {ta}) I 

Ui<j[U - tj) Uj<Mi - Xj) « \oti J 

Here the sign (— 1) M can be absorbed into the product Y[i=tj(K—tj) = (— 1) M Yli^jiti — Xj). The 
orthogonality of two different Bethe eigenstates was shown in ref.(||). From this expression 
taking the limit Aj — > ti one can easily obtain the formula for the norm of the Bethe eigenstate: 

where f(ti) = Yia^i c(t a — ti)/c(ti — t a ). Note that the matrix = —d/dtj\n(a(ti)/f(ti)) 
can also be represented in the form: 

2rj 

N* = ~hi(a(t i )))-T, 2 '' 




dU ^ (t ai + 7]) (t ai - rj) ' 

where tij = tj — tj. The first term in Nu should be understood as — (ln(a(t))) (t = ti). To 
calculate the correlators in the continuum limit it can be useful to represent the matrix iV as 
a product of two matrices as 



Nij = ]T (Sia + (1 - S^N^N^)- 1 ) (SajNjj) 



where Nu ~ R{ti) - the density of the variables ti in the continuum limit. The determi- 
nant expressions for the formfactors can be rather straightforwardly obtained from the above 
formulas using the explicit solution of quantum inverse scattering problem presented in ref. [|J. 

In conclusion, for the calculation of correlation functions using the formfactors, it seems 
necessary to develop the method to evaluate the dependence of the formfactors of local oper- 
ators {crf' z ) on the quantum numbers (quasiparticle numbers) characterizing the ground and 
the excited states, the problem which was not solved by now. The physical results have been 
obtained only for the simplest formfactors of the operator a z for the massles and massive 
|T2| (Baxter's formula for spontaneous magnetization) regimes of the XXZ chain. 
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We give in this appendix the proof that the operator (|13"D introduced in the text is the 
factorizing operator: 

6 = F 1 F 2 ... F N , Fi = (1 - h^ + T^i, (23) 

As it was explained in section 2 it is sufficient to prove the property O = Si+nO^' %+1 \ where 
by the definition 6^ i+1) = 6 12 .. .i+i,i,..N f° r the one particular permutation (i, i + 1). Consider 
the component FjFj+i in eq.fl23|). Substituting in this term Tj = Si + i^Ti l+2 '" N \ where 

rp(i+2...N) _ q q 
1 i — ^i+2,i ■ ■ ■ •JNii 

we obtain: 

FiF i+1 = (1 - hi)(l - h i+ i) + Si + i^ t+2 '-- N) T^ 1 2 ''' N) h i hi + i + 

S l+l ,Tt 2 - N) hi{\ - h l+1 ) + TlX 2 - N) h l+l {l - hi). (24) 

Moving Si + i t i to the right in the second term of eq.(^) using the Yang- Baxter equation, 
using the unitarity of the S- matrix and taking into account that S it i + ihih i+ i = hih i+ i we 
rewrite ( |24"D in the following form: 

Si+i.i ( (1 - hi)(l - h i+1 ) + S i ^ i+1 T l { ^ 1 2 - N) T l it+2 - N) hihi +1 + 

S iM1 Tl +2 - N) h l+1 (l - fa) + Tt 2 - N) hi{l ~ n i+ i) ) (25) 

Comparing eq.(|2^) and eq. (PB|) , we see that the expression in the brackets in ( P^| ) coincide 
with the right-hand side of (p^) up to the interchange of the indices % «-> i + 1 (together with 
the corresponding parameters £j, ^i+i). In the expression (^3j) the indices i, i + 1 are contained 
only on the left from the operator FiFi + \ in the operators T n for n < i: 

Thus moving the matrix Si + \^ in (|25| ) to the left and using the Yang - Baxter equation, we 
interchange the indices i, i + 1 in the whole expression for the operator O (p3|) and arrive 
exactly to the formula flUl) . 



Appendix B. 

In this appendix we derive the expression for the function $m(^, t) - the partition function 
of a six- vertex model with domain wall boundary conditions, and point out its properties 
used in the text. For simplicity we consider the rational case. We define this function using 
the operators B(t) corresponding to the monodromy matrix defined according to Tq(£) = 
Soi 5*02 • • • Sqm with the inhomogeneity parameters £i, . . . £m (that corresponds to its definition 
in the text): 

$m({£}, {t}) = ({lll..l}Af|B(fi)fl(* 2 ) . . . B(t M )\0). 
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It is just the same function that appears in the formula (p0|). The function is determined 
entirely by the following properties. 

1) At M=l, $i = 77/(t-f + »7). 

2) <&m({£}> W) i s a symmetric function of the variables {t} and separately {A}. 

3) If one chooses the initial (non -unitary) normalization of the S- matrix (and the cor- 
responding B operators), $m is a polynomial of degree M — 1 in each of the variables tj, Aj. 
The polynomial in is completely fixed by its values at M different points (for example at 
the points 

4) At t\ = £i the function reduces to $M-i({6i}i^i> an d the same for all the 
other pairs of variables. This can be easily seen using the symmetry property, by placing the 
corresponding variable £j to the end of the chain (at M-th site). Altenatively, one can use the 
representation of B- operators in the F- basis (|TE|). 

These properties determine the function <& M unambigously (see eq.(|21~D in the text). One 
should mention that it does not have poles at t{ = tj and & = £j since at these values the two 
lines or coloumns in the determinant coincide. 

One can also determine the function $m(£,£) directly using the operators in the F- basis. 
This was done in ref. ||. Let us present here this proof with some technical modifications which 
could be important for the other problems of this type and will be used in the Appendix C. 
Consider the determinant 



det(M ii ; 




A,)(/, - Xj + '/) 



First, we modify the first line of the matrix M\j adding a linear combination of the other lines 
so that the determinant det M is unchanged: 

M[ 3 = M : j + ]T CM ,,, 
where the coefficients are chosen to be 

c TT (tx - A/3 + Tj) -pj- (h ~ t a ) 

X \ lih-Xp + r)) J£ x {t x -t a y 

It is necessary to calculate the following sum: 

V C x M xj = - n^i(*i-^) y> f[t ) tt ( — \ — \ 1 1 (26) 

where 

/(*.) = n t'^'Vl 
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To calculate the sum in eq.(^) consider the integral in the complex plane over the circle of 
large radius which is equal to zero due to the behaviour of the integrand at infinity: 

dz . M 

[Z - h)(Z - Xj) UaM Z ~ 

where 

f(z) = rr {z -^ + r]) 

The sum of the residues at z = t a is equal to the sum in (^j), so it is easily calculated to be 



(il — Xj) n a ^l(il — ta) (Xj — ii) X[ a ^i(Xj — t a ) 

Thus combining all terms we obtain the following expression for the new elements of the first 
line: 

13 (*i " Xj) (ii - Xj + V ) U (Aj - O ^ (ii - A/3 + 77) ' 1 j 

The same formula could be obtained for the first column of the matrix My. Note also that 
taking the coefficiens C x in a different way we could obtained the other expressions of this 
type. For example, taking the coefficients 

a^l a l A l ~~ C «J /3^M _ A /3J 

we obtain the new expression of the first column of the matrix M 1 in the form: 

M , = 1 TT (j» -t a + T]) (Ai - Xp) 

a (t t -X 1 )(t t -X 1 + r ] )^ i (Ai-i a ) ^ (U - Xp + rjY 

which is different from (p7). Note also an interesting property of the function <&m(£, i), which 
can be seen from these formulas. Namely the function $^ is equal to zero at the points 
£i = £j ± 77. We will use this property in the Appendix C. 

Now the derivation of the formula (£HJ) is straightforward. Using the operators B(t) in the 
F- basis and acting on the vacuum state by one of the commuting operators, say the operator 
B F (ti) first, we get the formula 

=Eh(ti)mt), Mh) = bih-d) n L _ 7\ 

i=l ajii C v^« 

(we change the arguments in comparison with eq. (|16|) due to our definition of the function 
<&at)j where the function F(£, t) corresponds to the action of the other operators B F {tj). Each 
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term in this sum corresponds to the occupied site i, so due to the form of the operators B F 
fll6| ) the action of the other operators 



B F (t M )...B F (t 3 )B F (t 2 ) 

on this state does not contain the variable Thus we single out the dependence of the 
function on the variables & and ti, and the function F(£,t) is equal to 

F(£,t) = ^-i(feWi,W«^i), 

and does not contain these variables. Apart from the factor in front of the determinant that 
procedure corresponds to the first line development of the determinant in (pip . One can check 
that starting from the last formula we get exactly the first line development of the determinant 
with the matrix elements M' u instead of Mu. In fact, assuming that the function $>m-i is 
given by the formula fl2ip we obtain from the last formula: 

*m(u) = n (t U -^u% _^ E(-irdet(M a ,)iff-^(e,t) 

lli<j\ l i 1 11/ iVSJ i 
V TT (*1 ~ t a ) yr (& ~ & + V) 

(*i - eo(*i -ti+v) :k (6 - to) & - & + 7?) 

Comparing this formula with (^7|) we see that the last term in the brackets is exactly equal to 
M' u and we obtain the exactly the determinant of the matrix M' . Since we have shown that 
the two determinants are equal this procedure gives us the recurrence relation which proves 



the formula (ZI). 



Appendix C. 

In this appendix we will establish the direct correspondence between the calculations of 
the correlators in the F- basis and the previous calculations [§J based on the commutation 
relations (or, alternatively using the method of the present paper). Namely, starting from 



the expression of the scalar product in the F- basis we obtain the general formula (20) for an 
arbitrary parameters {£}, {A}. To obtain the expression of the scalar product in the F- basis 
we insert the complete set of states labeled by the coordinates {x} = (xi, . . as follows: 

Sm = E (0|C F (A!) . . .C F (A M )|{x})({z}|B F (ti) • • • B F (t M )\0). 

{x} 

One can use the expressions for B F , C F from Section 3 to find the expression for Sm (which 
is possible due to the quasilocal structure of these operators). However the simplest way is to 
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reduce the matrix elements to the scalar products of the usual form: 

S M ({X},{t})= E ( II np \ p A (0\C(X 1 )...C(X M )B(U---B(UM 

xi,..x M a^Xj C y^ai SXjJ J 

(0\C(Z xl )...C(Z xu )B(t 1 )...B(t M )\0). 

This formula is obtained by transforming the operators B and C back to the initial basis and 
using the formulas of Section 2. The last scalar products can be easily evaluated using eq. (|20|) . 
Proceeding in this way we obtain the following expression: 

M M 

SM({X},{t}) = Ua(t i )l[a(\ j ) 

i=l j=l 

£ n ^y (*" ( « j n^)) n^) . w 

where we have introduced the sets of parameters {£} = U with = (£ Xl . . . £ XM ) 
and the products are over all elements of the corresponding sets of parameters 
{£}, {A}). The function given by the product of two brackets in eq.(p8|) has only simple poles 
at £ Xi = tj and £ Xi = Xj. We remind that the function $m(A, £ x ) Yl(£ x — A + 77) has no poles in 
£ Xi . The behaviour of this function at the infinity is ~ for each £ Xi . To take the sum we 
apply the same method as in the Appendix B. The problem is that now we have to sum over 
M coordinates x±, . . . xm with the condition Xi 7^ Xj. First, as an example, let us consider the 
following simplified sum which is quiet similar to that in eq.( 



s = e I n a ~ir~\ n 1 



Xi,..xm Qi^Xj i,j i^ x i 

To calculate this sum let us multiply both the denominator and the numerator by Yli^j(£ Xi — 
£ Xj ). Then we obtain the following sum: 



XI, ..Xm 



Yiaj^xi (Ca Cxi) • • • Ila^xM (£ Q Cx M ) i j i^Xi Xj) 



where the sum is still over the configurations with Xi 7^ Xj. However one can relax this 
condition since in the denominator we have the function which is not equal to zero at X{ = Xj, 
while the numerator is equal to zero at Xi = Xj. Then the sum is represented by the product 
of the sums of the type fl2"6|) (with some degree of £ Xi in the numerator) which can be easily 
calculated by the same method. In fact, by means of the polynomial 

II(& _ = E Cki,...k M £,i 1 ■ ■ ■ Cm 1 , 

i^j k lt ...k M 
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(ki G Z + ) the sum is represented as 



s ~ E Ck U ...k M (^2 •))■■■ (E tit i i - 

ki,...k u \xi ilsxj/ \x M Jy<,x M )/ 

where 

M 

/(6*) = II (& C- ) H(A, -O ). 

Each of the sums over Xi can be rewritten as in the Appendix B: 

e t|h = e i^p-v /(v) = m» - k) n ( A i - A 



Then we obtain the product of the other sums of this type and the procedure can be repeated 
in the opposite direction. Thus we obtain the following simple result for this sum: S = 
Ha=i Ui=i(^,a — K)^ 1 - Now using the same method one can calculate directly the sum (|28|) . 
Schematically, for each of the variables £ Xi one can take either the pole at £ 4 or the pole at 
Xj. Taking the pole at ti we get l/a(i») which cancels the corresponding factor in front of the 
sum in (p8[). Using the notations of the formula (pO]), that means that the corresponding ti 
belongs to the set {/i}. We should also obtain the factor c _1 (/i — v) arriving to the formula 
). In fact, using the same method as before we represent this sum as: 

1 \ 




11(6* - life, - ^ + vr^M(t,Q^ M (\Q 

The function in the second line of this formula F(£ x ,t, A) - is the function of 3M variables. 
The function F has the following properties. 1) F is equal to zero at X, = Xj, so we can relax 
the condition Xi ^ Xj and sum over all Xi independently. 2) The function F(£ x ,t, A) does not 
have poles at the points £ Xi = £ x . ± rj since the function $m(£)£c) is equal to zero at these 
points. This can be seen from the representation given by the last equation of the Appendix 
B. The only poles contained in the function F are at £ Xi = tj — rj and £ Xi = Xj — rj. 3) The 
behaviour of F at \£ Xi \ — > oo is ~ We can take the sum consequently over x±, . . 

Taking the sum over x± we obtain residues at some £j or Aj (which corresponds to the set 
{/x}). One can see that the residues corresponding to the poles contained in the function F 
are equal to zero due to the terms of the first line of the last formula. Then we repeat this 
procedure for x 2 , . . . %. Note that the pole we take for x 2 should not coincide with the pole 
we took at the previous stage (for x±). One could formulate the same procedure in another, 
equivalent way. Consider the function 

F'(£ x , t, A) = t, A) n(£* -t + v) U(Zx - A + rj) = 
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n cfe, - txs)*M(t, q$ m (\ &) n(e^ - <j- + 77) nc^ - a,- + V ) 

This function does not have poles at all and therefore is a polynomial. Then the sum is 
represented as the product of independent sums (as for our simplified example). Thus, we 
finally obtain the result: 

M l l 
E n-7-Tr U , $M(t,/*)$Af(A,/z), 

m...n M i=l u \h L %) ij HA 4 * u j) 

where the notations are the same as in the equation (p0| ). Thus the equation (|20f) is reproduced. 
Note that one could start from the equation ( p0|) and derive the equation (p8|) using the same 
method substituting Jli = Ili Ilj a>(K)/ Hi and taking the sum over /ii, . . . 
with a(/x) = Ila=i c(£ a — /i) (taking the residues corresponding to the poles at the points 

= 

In conclusion, let us present two equivalent ways to rewrite the expression (^0|) for the case 
when the parameters {t} satisfy the Bethe Ansatz equations. Expressing the functions a(tk) 
and using the explicit form of the function <3>m we get: 

s M (x,t) = - 1 1 A , E(-i) Pfc (-i) p -nKAn)) 

det(4, \p) det(A n , t a ) Y[(X P - A n + rj)(t k -t a + rj)(t a - A n + 7]){\ p -t k + r]), 

where again the products is over all elements of the corresponding sets of parameters. Here 
the sign factors P k , P n depend on the sets of the coordinates k%, . . . k m , n 1; . . . n M _ m in order 
to obtain the factors in front of the sum from the terms n(^fc — t a ) and EI (A/3 — A n ). We also 
used the simplified notations for the determinants ([21]). One can further rewrite this formula 
as follows: 

5W (a, t) = 1 1 , E(-i) P H-i) p "(-i) Mm IlKAn)) 

IK* - An + V) U(t - A, - v) det(4, A,) det(A n , t Q ) R {r^A II i ^'^V ) ' 

where we denote for example Y\(t — X n + v) = Uj Hi=i{U ~ X nj + rj). The other notations are 
the same as before. In this formula the poles at ti = Xj are contained in the determinants. If 
the two last products are dropped out this formula becomes the result of the decomposition 
of the determinant in eq.(E2). Then one can see explicitly that the residues in the poles in 



the formulas (|20"D and (^) satisfy the same recurrence relations. 
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